Anomalous Hall effect for the phonon heat conductivity in a paramagnetic dielectric 
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The theory of anomalous Hall effect for the heat transfer in a paramagnetic dielectric, discovered 
experimentally in [lj, is developed. The appearance of the phonon heat flux normal to both the 
temperature gradient and the magnetic field is connected with the interaction of magnetic ions 
with the crystal field oscillations. In crystals with an arbitrary phonon spectrum this interaction 
creates the elliptical polarization of phonons. The kinetics related to phonon scattering induced by 
the spin-phonon interaction determines an origin of the off-diagonal phonon density matrix. The 
combination of the both factors is decisive for the phenomenon under consideration. 

PACS numbers: 66.70.+f, 72.15.Gd, 72.20.Pa 



A novel interesting phenomenon is found experimen- 
tally in recent paper The matter concerns an analog 
of the anomalous Hall effect (AHE) for the heat conduc- 
tivity of an ionic paramagnetic dielectric. In fact, ap- 
plying magnetic field B in the direction normal to heat 
flow j, the authors have discovered an appearance of the 
heat transfer in the direction normal to B and j. For the 
complete lack of free charged carriers and negligible role 
of the spin-spin coupling at the parameters concerned, 
the transverse flow is naturally associated with the evo- 
lution of the phonon system. However, the magnetic field 
does not directly act on phonons and only polarizes para- 
magnetic ions. It is the coupling between phonons and 
subsystem of isolated ions carrying magnetic moment M 
that determines the formation of the picture observed. 
Note in this aspect that the phenomenon concerned is 
an analog of the AHE in the paramagnetic phase of a 
ferromagnetic above the Curie point, discovered a long 
time ago [1| ■ An existence of the effect found in [l[ has 
been confirmed in work '3] . In both d and i the choice 
of terbium gallium garnet (TbGG) was not accidental. 
In this compound for T < 10K the heat conductivity 
due to spin-phonon coupling (SPC) proves to be about 
two orders of the magnitude smaller as compared with 
gadolinium garnet 0] of the identical structure. 

As it is known, in most of ion dielectrics the spin- 
phonon relaxation is determined by the two-phonon pro- 
cesses, see, e.g., [1 and @. Usually this supposes the 
presence of the Kramers doublet in the ground state and 
the transitions via virtual excitation of higher Kramers 
doublets. Trivalent ion Tb 3+ has the even number of f- 
electrons and the Kramers degeneration is absent. How- 
ever, for rare earth ions with the even number of f- 
electrons, an appearance of the quasi-doublet structure is 
typical for splitting the multiplet in the crystalline field 
with the level spacing £12 of about several Kelvins in 
the ground state. TbGG is a striking example of the 
similar compound, see, e.g., Q. From the gene ral sym- 
metry reasons and the direct analysis one can 
conclude that the SPC under conditions concerned re- 
duces to a scalar production of magnetic moment M and 



the orbital moment related to the motion of surrounding 
ions. At low temperatures when the long wave acoustic 
branches alone are excited, the displacement of all atoms 
in the elementary cell for each mode is the same. One can 
show in this case that the total orbital moment of atoms 
in the elementary cell reduces to the angular moment of 
the cell center of gravity. We will assume that the mag- 
netic ordering temperature T c <C T and the magnitude 
of magnetic field allows us to neglect the dipole-dipole 
interaction between ions. This gives a possibility to an- 
alyze the phenomenon adequately, considering a simple 
model system with one magnetic atom in the elemen- 
tary cell and assuming high symmetry of SPC. Thus, the 
Hamiltonian of SPC can be represented as the spin-orbit 
interaction (see, e.g., [8() 



Hi =g^2,s n {U n x P„), 



(1) 



where U n and P n are the vectors of displacement and 
momentum of the center of gravity in n-th elementary 
cell and s n is the isospin typical for the quasi-doublet. 
For the isotropic SPC, one has s n \\ M n . Assuming that 
the spacing from the ground quasi-doublet to higher lev- 
els is As 3> T, magnetic moment M n is determined by 
the parameters of the quasi-doublet components with the 
regard of their occupation. In the approximation linear 
in Hi we can replace M n and, correspondingly s„, with 
the quantities M and s averaged over the crystal. 

The present paper is devoted to determining the heat 



conduction tensor 



odd (j^) antisymmetric and odd m 

magnetic field B [11]. The heat conduction problem is 
solved at low temperatures when the heat is carried with 
long wave acoustic phonons. A fraction of heat flow di- 



rected along 



M x VT 



arises from the elliptic renormal- 

ization of the polarization vectors of acoustic phonons 
with involving ([T]). The renormalization of the sound 
velocity does not take place. The effect occurs at an ar- 
bitrary character of the phonon dispersion spectrum and, 
in fact, at an arbitrary mechanism of phonon scattering 
and phonon-phonon coupling. In this communication we 



2 



restrict ourselves with considering phonon scattering at 
magnetic ions with the involvement of the level system 
which appears in the crystal field. The elements of den- 
sity matrix p ss ,(k) = \a ks a ks ') off-diagonal in branches 
and found from the generalized Boltzmann kinetic equa- 
tion play a key role in the formation of x°^ dd - 

In recent paper [13] there has been attempted to ex- 
plain theoretically the phenomenon observed in [l[ . Tak- 
ing SPC in the form analogous to |T]), the authors have 
chosen a simplified model for the phonon spectrum in 
which two coincident transverse modes always exist for 
an arbitrary wave vector. Such degeneration is decisive 
for the consideration on the whole. In fact, in the general 
case the acoustic branches are not degenerated in most 
part of phase volume. The regions in which one can meet 
degeneration occupy, as a rule, relatively small fraction 
of the phase volume. This means a necessity of the ki- 
netic consideration with involving the phonon scattering 
mechanisms. 

2. Treating the phonon system in the harmonic 
approximation, we write the general Hamiltonian of the 
system with (p} as 

H = H +H 1 , Ho = Y,^Pn- 7J YJ2 D nn'KU b nl . 

n 7in' 

(2) 

Here mo is the total mass of the elementary cell. Let us 
determine the Hamilton's equations from ([2|). For veloc- 
ity V®, one has straightforwardly 



limit. The latter ones can be chosen as real. In the ap- 
proximation linear in the SOI one can conclude from (5) 
that the phonon spectrum remains unchanged and the 
polarization vectors alone renormalize. Let us introduce 
notation 



U^(k) = e^ a (k) + 5e a s (k). 



(7) 



Then from equation ([5]) one finds the small rotation of 
the polarization vector for s-th branch 



Sei 



„. r ((e 2 xjQj) 
-2igwi[- r ^ — e 2 



((eg x ei) s) 



-e 3 j 



(8) 



Accordingly, for s = 2 and 3. In the linear approxima- 
tion in SOI the condition of orthonormality holds for as 

p a *e a , = S , 

To find the phonon energy flow j c with the presence of 
SOI, one can employ the general results obtained in [l3j ]. 
The successive use of the continuity equation, reflecting 
the local law of energy conservation for the slow spatial 
temperature variation, results in the expression 



J 



2V 



m 



Er>c jjab jjayb 
^nn' nn' ^ n * n' ' 



(9) 



Hereafter the volume is V=l. In this expression the 
velocity V£, is determined by relation © . Let us expand 
vectors £7™ (and V% — U%) into the normal modes in the 
representation of secondary quantization (% = 1) 



V: = U a n =P«/m +ge abc s b U c n , 



(3) 



where e a bc is the antisymmetric third-rank unit tensor. 
Using explicit expression ([3]), the equation of motion can 
be transformed to the form 

-K=J2 Dt> K - U b n >) + 2. 9 e abc l>y . (4) 



Assuming further the consideration in the linear ap- 
proximation for the spin-orbit interaction (SOI), we re- 
placed P b with molI b in the last term. A usual Fourier- 
transformation of linear equation ([4]) results in the dis- 
persion equation 



where 



D t = H D i ( J - e ^ E D t = °- 



(5) 



(6) 



Here R = r n — r n i , r n is a position of the center of gravity 
in 7i-th cell. In zeroth order in SOI the solution of a set 
of equations (|5]) determines the dispersion law oj s (k) for 
three acoustic branches and, correspondingly, three or- 
thonormal polarization vectors ei°^ a (£;) in the long wave 



E 

ks 



i 



2m Nu k . 



■ exp 



(ikr n ) [ e -{k)a ks +eT{-k)a + _ ks 

(10) 

Substituting these expressions into © and employing re- 
lation ^ , we find for the averaged operator of the energy 
flow 

4 £— ; V uj ks > V Lu ks 

(11) 

where p ss ,(k) — (ci ks O'ks'/- Keeping the approximation 
linear in SOI, we neglect anomalous averages {a- ks a ks >) 
and (a ks a- ks i) in (|lip . For zero order approximation 
in SOI the polarization vectors in (fl"Tj) are real and the 
expression in the figure brackets is symmetric over in- 
dices s,s'. Thus only the symmetric and real component 
of density matrix gives nonzero contribution to ((jlip). 
Assuming p ss ,{k) — n s (k)5 ss ' and the relation 



{VlDf ) elAs = dA ks /dk c = 2u ks C c ks 



(12) 



where C s — duj s /dk is the sound velocity, from (jlip we 
find an ordinary expression for the heat flow 



O'o) = ^2^ksC c ks n s {k) 



(13) 



ks 
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The contribution linear in SOI appears in (fTTj) with 
involving the renormalization of the polarization vectors 



a* b 
e ks e ks 



+ e 



(0)o 



Set 



(14) 



Substituting relations {8| into this expression, we find 
that the contribution diagonal in modes vanishes. The 
off-diagonal term with ks — 1 and ks' = 2 equals 



OA„ ( g 2 xei)s / a „6 , , -a-f^ 



9 - ifsxealf 



w 2 e 1 e 3 + lig^i_J^LOie 3 e 2 



(15) 



With substitution lfl"4] ) .(fT5" J) into (fTTj) the expression m 
the figure brackets, proves to be antisymmetric over in- 
dices 1 and 2. The contribution linear in SOI in (fTTj) does 
not vanish if the off-diagonal terms of density matrix are 
nonzero. In the expression (|15p the first term plays a key 
role. To avoid the cumbersome expressions, we will find 
the approximate value for the flow (fTTj) , keeping the first 
term alone (j% Q ) = (j£ 2 ) + (jf 3 ) + (j^), where 



0*12 ) — i-9j2k (y^ + yW, 

2u 1 u 2 (e 2 'Xe 1 )s (/-to, r<a\ „ 
(w*- u *) \ U 2 ~ °1 ) 9X1- 



(16) 



To find the flow finally, it is necessary to determine 
the imaginary part of nonequilibrium off-diagonal den- 
sity matrix p 12 . We employ a conventional procedure to 
derive the kinetic equation, see, e.g., [14[. The general 
equation of evolution reads 



id tPkk' = {[ata k >,H Q +H'}) 



(17) 



Here H' is the Hamiltonian for the phonon scattering and 
Hq is the starting Hamiltonian of noninteracting phonons 
(f2j). To simplify notations, we introduce a generalized 
index for the mode k = k, s. In the frequency represen- 
tation 



Considering the region of low temperatures T ^ 6K, 
we assume that the kinetics is determined by the scatter- 
ing of phonons at this quasi-doublet. The higher levels lie 
at £j > 50K and we neglect rescattering through them. 
In this case Hamiltonian H' can be represented in the 
general form 



Tjl — __ \ "* \ " J{g-f)r n Aotafa + 

n - N 2^2^ e A fgt,na f ag. 



(20) 



Here a — 1,2 numerates the doublet levels and = 1, 0. 
The commutator in (fT5j) equals 

[afcOfc'^nO/ a fl]) = { a kda g )5k's - <^s 

Using equations (fT5j) for these correlators, we obtain 
some quadruple operator correlators. 

Restricting with the Born approximation, we decouple 
these correlators in the mean-field approximation. As a 
result, we have at uj — > 



Ikk' 



J_ V 1 p fl — P \ Aaa Aaa 
N l^qa^ay 1 r <* J^k' g^gk 



{5 (cjfe - w g ) {n k - rig) +5 {uj k > - u> g ) (n k > 



(21) 

where P a = is the occupation value of level a. From 
(f2Tj) we can conclude that I kk i = I k 'k- As a result, we see 
that off-diagonal matrix p kk , (I19p is antisymmetric. For 
the equilibrium distribution, l k y = 0. Nonzero result ap- 
pears only for nonequilibrium distribution f k — n k —N^ 
due to temperature gradient. In the r-approximation the 
expression (|21ll can be represented as 



(22) 



Ikk' — — i^kk'fk + ^k'kfk') , 



L0 



(23) 

The usual solution of the problem for the longitudinal 
heat conduction yields 



/ \ t t /r + zjh\ „o^ /fe = ~ (cJk/flkkT 2 ) Ni 0) fl + N^) C k VT (24) 

(w + w fc -uj k ,) Pkk , = I kk >, hk' = {[ala k r,H']) (18) ^ fc/ ** ' k \ ^ k J y ' 



The formal stationary solution (uj k ^ LJ k i) reads 



P kk ' 



kk ' 



Wfc — U)k> 



(19) 



At low temperatures, as we have already mentioned, 
the main mechanism of scattering is connected with the 
coupling between phonons and magnetic ions. The split- 
ting of rare-earth ion multiplet with the crystalline field 
V cr and the reduction of the symmetry for oscillating part 
V cr make decisive the Raman two-phonon scattering via 
virtual excitation of the other multiplet levels [H, Q . For 
garnets like Tb 3+ , the ground state is a quasi-doublet 
with the small spacing of about several Kelvins between 
the doublet components. 



The substitution of the off-diagonal density matrix 
(fTJjl into (fUj) with involving (|2"2" J) -([2"l j) solves the prob- 
lem for the heat transfer due to SOI. The relaxation 
frequencies (|23p depend on a product of the transition 
amplitudes with k ^ kl instead of a usual square of the 
modulus of the transition amplitude when everything is 
determined by the diagonal density matrix. This typi- 
cal feature is inherent in any scattering mechanism, in 
particular, due to phonon anharmonicity. 

Let us return to the notations adopted in (fTTjj) and 
make a significant remark about the symmetry for a 
product of polarization vectors with the same wave vec- 
tor k entering this expression. In the long wave approx- 
imation for the crystals of the sufficiently high symme- 
try the off-diagonal elements of the dynamical matrix 
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read Df 



k a k b , e.g., [15J. Rewrite the dispersion equa- 
tion ([5]) for zero order approximation in SOI in the form 



One can show that, 
with varying the sign of the separate projection of wave 
vector, D^ a and u> 2 s remain unchanged. Then from the 
given equation it follows that the projection of polariza- 
tion vector e£ s changes its sign for the substitution k a 
with (— k a ) conserving the sign of the other projections. 
This property can be expressed in terms 

e a ks = % (*) signk" (25) 

Here e" (jtj is a unit vector remaining unchanged with 
varying the sign of an arbitrary projection of wave vector 

Let magnetic field and magnetic moment M be in the 
direction of the z-axis, while the longitudinal VT be in 
the direction of the cc-axis. Then, substituting [)19p and 
(|22 p -(|25 p into we find, restoring H, k B and volume 
V, for the transverse component of the heat conduction 



"12 



l 23 



*31 ) 



"12 



B 9 — 



Efc sZ (ei x e 2 y 



2uJiU-2signk x signk v / y 



{u>\— w|)(o)i— w 2 ) 
F\2 = 



U> 2 ) 



(cf-cl)(F 12 cT 




(26) 



r(0) 



Hii (fe B T) 2 



(l + ivf) 



Recall that = (s z ). The expression obtained 

demonstrates an existence of the anomalous phonon Hall 
effect under comparatively general conditions conserving 
the intrinsic symmetry. It is interesting that the kinetics, 
reflecting phonon scattering character, enters (|26[) only as 
a ratio like O12/O11. This holds for any dominant mecha- 
nism of scattering. Under these conditions the transverse 
heat transfer is determined mainly by the spectral prop- 
erties of the phonon system. Note, as it follows from (f26]) . 
that a quantitative enhancement of the effect takes place 
if the region in which branches uj ks come close is notice- 
able. We will give a quantitative estimate of the effect, 
assuming that the role of these regions in the integral 
determining (|26p is limited in the general case. 

For the estimate, one should determine the magnitude 
of parameter g introduced formally. The SOI ([1]) due 
to two-phonon processes accompanying vibrations of the 
crystalline field has the same origin as Hamiltonian H' 
(|20p . The Hamiltonian of singlc-phonon interaction is 

written usually for a mode {k-^j m the simplified form 

V (Uks/a) (ka), where V is close to the magnitude of the 
static crystalline field (see, e.g., [6J]). Factor (ka) appears 
with regard to the relative motion of magnetic ion and 
surrounding atoms. In second order in this interaction, 
the standard unitary transformation results in the spin- 
orbit interaction as (see [H) 



(s z )V' 2 U ks V ks 



'ks 



To simplify, we introduce here an averaged value for 
sound velocity C and take uj ks ~ Lo ks > ■ Choosing two- 
mode contribution to SOI, written in the form {]]), we 
have (s z ) gmoUksVks' ■ In fact, we should mean the dis- 
placement in Hi as a relative displacement of ions. Thus 
in this expression, as well as in the final one (|26p . 
g = g (ka) should be implied instead of g. Comparing 
the expressions presented and assuming the temperature 
T > £12, we find g ~ V' 2 / (hm Q C 2 ). 

Returning to (|26p . let us make an integral dimension- 



less. Dividing sound velocity by C, we obtain ~ 
3tt- 2 (HC/T) J dkk 2 ... The remaining part of the inte- 
grand is dimensionless. Putting that (J- e ffB < ksT, we 
find within the accuracy of a numerical factor 

x«" ~ k B g (k B Q D /hC) (T/Q D ) (n eff B/k B T) . (27) 

Taking V ~ 5QK, C ~ 4 * 10 5 cm/s, D ~ 400^, 
we arrive at ~ 10" 7 (fi eff B/k B T) (W/cmK). Using 
the experimental value for [1,4] at H e ffB jk B T ~ 1, 
we find for the Hall angle 



77 = x yx /x xx ~ 2 * 10" 



(28) 



In reality, this ratio may be larger since the numerical 
coefficient omitted in ([28]) can take the value significantly 
larger than unity due to the character in the behavior of 
integrand (j2"6"j) . 

So, in the general case the anomalous phonon Hall ef- 
fect originates from a combination of two important fac- 
tors. The first of them is associated with an appearance 
of elliptic renormalization for phonon modes as a result of 
the spin-orbit interaction between phonons and paramag- 
netic ions. The second factor is connected with the proof 
that the off-diagonal density matrix of phonons does not 
vanish at the kinetics of the longitudinal heat transfer 
with the presence of temperature gradient VT. The lat- 
ter predetermines a necessity to solve nonconventional 
kinetic equation in the problem of phonon heat trans- 
fer. The model of a crystal and interaction used in the 
paper is undoubtedly simplified. That is why, the fact 
that the estimate obtained for the Hall angle ([28|) with 
using the real parameters proves to be comparable with 
the experimental magnitudes [l|, Q lying within interval 
seems very encouraging. 
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